Spin resonance peak in Na^CoC^ • 2/H 2 superconductors: A probe of the pairing 

symmetry 
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Motivated by the recent discovery of superconductivity in 
layered Na^CoC^ • 2/H2O compound, we investigate theoret- 
ically the spin dynamics in the superconducting state with 
three possible p x + ipy, d + id'-, and /-wave pairing sym- 
metries in two-dimensional triangular lattices. We find that 
a spin resonance peak, which is elaborated to have a close 
relevance to the relative phase of the gap function and the 
geometry of the Fermi surface, appears in both in-plane and 
out-of-plane components of the spin susceptibility \ for the 
spin-singlet d + id'-wave pairing, while only in the out-of- 
plane(in-plane) component of x f° r the spin-triplet p x + ip y - 
wave (/-wave) pairing. Because the pairing symmetry in this 
compound is still hotly debated at present, this distinct fea- 
ture may be used to probe or determine it unambiguously 
with future neutron scattering experiments. 

PACS number: 74.90.+n,74.20.Rp,74.25.Ha 



Recently, the superconductivity with T c ~ 5K is found 
in the CoCVlayered material, Na^CoC^ • 2/H 2 [1]. This 
compound consists of two-dimensional(2D) C0O2 layers 
where Co atoms form a 2D triangular lattice and is sep- 
arated by a thick insulating layer of Na + ions and H2O 
molecules. Although its T c is relatively low, it has still 
attracted much attention because it shares some similar- 
ities with high-T c cuprates, in particular it may be an- 
other unconventional superconductor evolving from the 
doped Mott insulator. Therefore, the investigation of 
this compound is expected to give insight on the mecha- 
nism of unconventional superconductivity. So far, several 
quite different theoretical proposals [2-7] for its pairing 
symmetry, such as the spin-singlet d + zrf'-wave and the 
spin-triplet p x + ip y -(or even /-)wave, have been put 
forward. On the other hand, experimental results on 
the pairing symmetry reported by different group have 
also been controversial, even with the same experimental 
method [8-11]. 

In unconventional d- and p-wave superconductors, the 
gap function A(k) changes the sign(phase) around the 
Fermi surface and thus would lead to zeros(nodes) in the 
superconducting(SC) energy gap |A(k)|. Therefore, one 
can in principle determine the gap function and conse- 
quently the pairing symmetry, by measuring the distri- 
bution of the phase and/or node positions. In practice, it 
is the node positions rather than the phase that can be in- 
ferred in usual thermodynamic, transport, and NMR ex- 



periments. Therefore, the probe of the node position has 
been mostly used in the clarification of the pairing sym- 
metry in unconventional superconductors [12]. However, 
the much debated pairing symmetries so far proposed for 
the SC state of Na^OoC^ • 2/H2O are the broken-time- 
reversal-symmetry d + id'- and p x + ip y - waves. In this 
case, the energy gap |A(k)| = ^/A^(k) + A 2 ,, (k) is node- 
less. So, a probe which is directly related to the phase is 
of special importance for the determination of the pair- 
ing symmetry in this compound as well as in other un- 
conventional superconductors. This is evident by noting 
the crucial impact of the tricrystal phase-sensitive experi- 
ment [13] on the determination of the dominating d x 2_ y 2 
pairing symmetry in high-T c superconductors. In this 
Letter, by noting that a spin resonance peak appears in 
different components of the dynamical spin susceptibil- 
ity x for all three possible pairing symmetries in the 2D 
triangular lattice with the nearest-neighbor(n.n.) pair- 
ing interaction, we show that the identification of the 
spin resonance peak in the SC state, which can be car- 
ried out by neutron scattering experiments, may also pro- 
vide an unambiguous clue to probe/determine the pair- 
ing symmetry in this compound. We elaborate that the 
occurrence of the spin resonance peak in a specific com- 
ponent of the dynamical spin susceptibility exclusively 
corresponds to a definite change of the phase of A(k). 

To handle both the spin-singlet and spin-triplet super- 
conductivity in the same model as well as to capture the 
essential physics of electron correlation, we employ a phe- 
nomenological t—U — V model [14,15] on a 2D triangular 
lattice, in which an effective n.n. pairing interaction (V) 
is responsible for superconductivity and an on-site Hub- 
bard U for the electron correlation. After choosing the 
mean-field parameter Afj = V(< c^Cji > ± < c^Cj-j > 
)/2, we can write the effective Hamiltonian as, 



<ij>,<J i 



(i) 



where the upper sign is for the spin-triplet pairing state 
and the lower sign for the spin-singlet pairing state. Note 
that the actual superconducting pairing mechanism is 
still unclear at present and thus is not specified in this 
phenomenological model. 
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In the triangular lattice, the dispersion of quasiparti- 
cles is 



e k = — 2t [cos k x + 2 cos — cos — -^-\ 



fi 



(2) 



For 2D triangular lattices with the n.n. supercon- 
ducting pairing interaction, only d x 2_ y 2 ± id xy -wave, 
p x ± ipy-wave, and /-wave pairing states may exist [16] 
: i) A d k +ld ' = A {[cos(fc x ) - cos(k x /2) cos(V3fc y /2)] + 
iV3sin(fc :c /2)sin(V3fcj / /2)}, ii) Af +lPy = A {[sin(fc x ) + 
sin(fc x /2) cos(v / 3fcy/2] -M\/3cos(fc x /2) sin(\/3fcy/2)}, and 
iii) A{ = A {sin(fc x ) — 2 sin(fc x /2) cos(v / 3fcj / /2)}. Given 
the attractive interaction V and with i < [2,4,15], the 
mean- field calculation of Eq.(l) shows that nearly de- 
generate singlet d + id'- and triplet /-wave solutions are 
favored at doping n = 0.4, while the triplet p x + ip y - 
wave is stable at n — 1.35, where n is the average elec- 
tron number per site [15,17]. In order to compare the 
results for different pairing symmetries with roughly the 
same SC gap, we have chosen V = 0.75t for n = 0.4 
and V = 1.7* for n = 1.35, which gives A = 0.015t 
for /-wave, A = 0.014i for d + id-wave (n = 0.4) and 
A = 0.15t for p + ip-wave (n = 1.35). The effective on- 
site Hubbard interaction is assumed to be U = 3t [18]. 
The bare spin susceptibility is given by, 



c £(fc.g)(*£ g -i) 

w - Q£„ + iT 



where the coherence factors are, 



C? j (k,q) = [l± 



e k e k+q + Rc(A k A* k+q ) 

Ek+qEk 



(4) 



for the spin-singlet pairing and if ij = zz (the out-of-plane 
component of x) for the spin-triplet pairing, and 



cg(fc, g ) = [i± CfcCfc+q ; Re( ^ ) ] 



(5) 



if = H — (the in-plane component of \) for the spin- 
triplet pairing. F± q = f{E k+q ) ± /(S fe ) and = 

E k ± with S fc = y/ei + |A fc | 2 and /(S fc ) the Fermi 

distribution function. Near T = 0, only the first term 
in Eq.(3) with the coherence factor C~, involving the 
creation of quasiparticle pairs, contributes to the spin 
susceptibility. An essential difference in the coherence 
factors between the spin-singlet pairing(or the compo- 
nent Xzz for spin-triplet pairing) and the component x+- 
for spin-triplet pairing is a sign difference in front of 
Re[AfcA£ + ]. This sign difference is the key point for 
our later discussions. 



The many-body correction to the spin susceptibility is 
included by the random phase approximation [18]. In 
this way, the renormalized spin susceptibility is given by, 



Xv (q, w) = x% (q, [i - Ux% (q, w)] " 



(6) 



The momentum dependences of the dynamical spin 
susceptibility for w = 0.02t in the SC state (T = 0.0001*) 
are presented in Fig.l. For the p x + ipy-w&ve pairing, 
a peak near Qi = 0.94 x (2tt/3,2tt/\/3) can be seen. 
While, for the /-wave and d + id' -wave pairings, a peak 
near Q2 = (0, \/37r/2) appears. An obvious feature, seen 
from the figure, is that these peaks depend only on the 
doping density and is irrespective of the symmetry of the 
pairing spin state and the components of Imx- This is 
due to the fact that these peaks arise from the nesting 
of Fermi surface which is determined only by the doping. 
The nesting wavevectors Qi and Q2 which correspond 
to these two peaks are denoted as dashed lines with ar- 
rows in Fig. 2. We note that a much sharper peak occurs 

around q = (0, 0) for the in-plane component Imxn in 

the case of the spin-triplet pairing. This peak already 
presents in the normal state (not shown here) and reflects 
an enhanced ferromagnetic fluctuations which may arise 
from the substantial density of state at the Fermi level. 
However, it is suppressed highly for spin-singlet pairing 
and for the out-of-plane component in the spin-triplet 
pairing as shown in Fig.l. 

In Fig. 3 (a), we present the frequency dependence of 
Im% at Qi for the spin-triplet p x + ip y -wave pairing. It is 
seen that a spin resonance peak occurs near uo — 0.05t for 
the out-of-plane component Imxzz, but is absent for the 

in-plane component Im%^ However, in sharp contrast, 

the spin resonance peak appears in the in-plane compo- 
nent of Imxn rather than in the out-of-plane component 

for the spin-triplet /-wave pairing(Fig.3(b)). Moreover, 
the spin resonance peak can be found in both components 
for d + id' pairing, via the relation Imx+- = 2Imx zz 
which holds for a spin-singlet state. These distinctly dif- 
ferent features for all three possible pairing states are sig- 
nificantly important because they may be used as an un- 
ambiguous clue to probe/determine experimentally the 
pairing symmetry in this compound. 

To understand our observation clearly, we plot the 
bare spin susceptibility Im^o m Fig. 4. It is clear that 
a peak is evident at the spin gap edge followed by a step- 
like decrease just below the gap edge in the channels in 
which there is a spin resonance peak. Using the Kramers- 
Kroenig relation, we will obtain a logarithmic singularity 
in its real part Rexo- Thus, the RPA correction will fur- 
ther magnify this effect and lead to a sharp peak near 
the gap edge. This indicates that a peak just above the 
spin gap edge is the source of the spin resonance. 

From the BCS theory, the density of statcs(DOS) is 
divergent just above the SC gap edge, which is expected 
to affect some physical quantities. However, the effect is 
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limited by the coherence factor. The most evident effect 
of the coherence factor [Eqs.(4) and (5)] is for energies Ek 
and Ek+ q near the quasiparticle gap edge A , in which 
it is either ~ or ~ 1, depending on the relative sign 
of Ak and Ak+ q - Specifically, C~ is negligible unless 
Afe and Ak+ q are of opposite signs for the spin-singlet 
pairing [19] and for the out-of-plane component in the 
spin-triplet pairing, or of the same sign for the in-plane 
component in the spin-triplet pairing. With these general 
considerations, let us now address why a peak appears in 
one component of the spin susceptibility, but is absent in 
the other. In Fig. 2, we plot the phase(+ sign denotes the 
phase 0, — sign the phase it) and node position(dotted 
lines) for various terms of the three possible pairing sym- 
metries. The Fermi surface for n — 1.35, where the 
Px + «Pj,-wave is favored, is a circle around (0, 0) point. 
For cither the p x or p y term, the two half circles sepa- 
rated by the line nodes will have the opposite(different) 
sign(phasc) of the gap function A&. Therefore, for the 
wave vector Qi = 0.94 x (27r/3, 27r/\/3), A& and Ak + Q 1 
have the opposite sign. According to Eqs.(4) and (5), 
the coherence factor C~ is appreciable for ij = zz and 
vanishes for ij = + — . As a result, the DOS peak shows 
up in the out-of-plane component ImXzz and does not 

in the in-plane component Im%_| , as shown in Fig. 4(a). 

Also for the spin-triplet pairing but with the /-wave, A& 
and Afc + g 2 connected by Q2 = (0, y/3n/2) are of the 
same sign(phase)[Fig.2]. Therefore, the DOS peak ex- 
ists in Imx^ , instead of in Im^zz [Fig. 4(b)]. The most 

definite demonstration of this argument can be found in 
the case of the d + id' wave, where an appreciable coher- 
ence factor C~[Eq.(4)] requires that A^ and Aj-+ q have 
the opposite sign. However, from Fig. 2 one can see that, 
though the A^'s connected by wave vector Qi satisfy 
the requirement for the d' term, those for the d term 
do not. To see their effect, we have calculated the re- 
sults for d and d' terms, separately. As shown in the 
inset of Fig. 3(b), we find no peak for the d-wave term, 
but a sharp peak for the d'-wave term. Remembering 
that the term Rc[A^ d ' A&f *] = A^ +g + Af A{ +q , 
one will expect that it is the effect of d' term dominates 
for the d + id' wave. The only difference between the 
cases of d- and d'-wave pairings is their sign (phase) of 
the gap function. So, it shows definitely that the spin 
resonance peak depends uniquely on the relative phase 
of the gap functions connected by the transition wave 
vector, when the Fermi surface is given. Therefore, the 
spin resonance peak appears in specific components of 
the dynamical spin susceptibility with distinctly differ- 
ent ways for the proposed possible pairing states. We 
suggest that neutron scattering experiments can probe 
the above mentioned spin resonance peak, as done for 
high-T c cuprates [20], and thus to identify the pairing 
symmetry of Na K Co02 • yH 2 superconductors. 

Before concluding the paper, let us use the above idea 



to address the anisotropic suppression of the spin re- 
sponse at q = (0,0) shown in Fig.l. At q ~ (0,0), the 
two gap functions connected by q will surely have the 
same phase. So, the coherence factor C~ will be neg- 
ligible for the spin-singlet pairing and the out-of-plane 
component of the spin-triplet pairing, but it is not for 
the in-plane component. Thus, comparing with that in 
the normal state and in the latter case, the spin response 
in the former case is strongly suppressed. 

In conclusion, we have found that the spin resonance 
peak exists in distinctly different ways for three possible 
pairing symmetries proposed for newly discovered C0O2 
layer materials Naa;Co02 ■ yEbO, and suggested to use it 
as an unambiguous clue to probe/determine the pairing 
symmetry in future inelastic neutron scattering experi- 
ments. Moreover, we have elaborated that the spin reso- 
nance peak has a close relevance to the relative phase of 
the gap function and the geometry of the Fermi surface. 

After the completion of this work, we noticed a 
preprint [21] by Li and Jiang, in which a spin resonance 
peak was predicted only for the spin-singlet d + id'-wave 
pairing in the framework of the t — J model on triangular 
lattices and by using the slave-boson approach. 
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and the RGC grant of Hong Kong. 
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FIGURE CAPTIONS 

Fig.l Momentum dependence of the spin susceptibility 
Im% with to — 0.02t for (a) p x + ip y -w&ve at n = 1.35, 
and (b) /- and d+id'- waves at n = 0.4. The momentum 
is scanned along the path shown in the inset of Fig. 1(a). 

Fig. 2 Fermi surface (think lines) and phase (±) of the 
gap functions for three possible pairing symmetries. The 
dotted lines denote the node positions which separate the 
regions with different phases, and the dashed lines with 
arrows represent the transition wave vectors Qi and Q 2 
as indicated in Fig.l. 

Fig. 3 Frequency dependence of the spin susceptibility 
Imx in the SC state (T = O.OOOlt): (a) for the p x + ip y - 
wave at n = 1.35 and Qi, (b) for the /- and d + id'- 
waves at n = 0.4 and Q 2 . The inset of Fig.(b) shows the 
results for the pure d and d! waves, respectively. 

Fig. 4 Frequency dependence of the bare spin suscepti- 
bility xo in the SC state (T = O.OOOlt) calculated from 
Eq.(3): (a) for the p x + ip y -w&vc at n = 1.35 and Qi, 
where the solid line indicates the out-of-plane component 
and the dashed line the in-plane component x+_; (b) 
for /- and d + id'- waves at n = 0.4 and Q 2 , where the 
solid line indicates x° 2 of the / wave, the dashed line x+- 
of the / wave and the dotted line x° f° r d + id' wave. 
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